Effects of impurities on superconductivity in noncentrosymmetric compounds 
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We microscopically derive the Ginzburg-Landau free energy functional for a noncentrosymmetric 
superconductor with a large spin-orbit splitting of the electron bands, in the presence of nonmagnetic 
impurities. The critical temperature is found to be suppressed by disorder, both for conventional and 
unconventional pairing, in the latter case according to the universal Abrikosov-Gor'kov function. 
The impurity effect on the upper critical field turns out to be non-universal, determined by the 
pairing symmetry and the band structure. In a BCS-like model, T c is not affected, while H C 2 
increases with disorder. For unconventional pairing, both T c and H C 2 are suppressed by disorder. 
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I. INTRODUCTION 

The recent discovery of superconductivity in a heavy- 
fermion compound CePt 3 Si, see Ref. [l|, has renewed in- 
terest, both experimental and theoretical, in the prop- 
erties of superconductors without inversion symmetry. 
The list of such materials has been steadily growing 
in recent months and now also includes Ulr (Ref. (3), 
CeRhSi 3 (Ref. |), CelrSis (Ref. !), Y 2 C 3 (Ref. i), and 
Li 2 (Pd 1 _ X) Pt !C )3B (Ref. i). 

A peculiar property of noncentrosymmetric crystals is 
that the spin-orbit (SO) coupling qualitatively changes 
the nature of single-electron states, namely it leads to the 
lifting of spin degeneracy and the splitting of the energy 
bands. This has important consequences for supercon- 
ductivity. If the typical band splitting Ego is smaller 
than the superconducting critical temperature, then the 
effects of the SO coupling can be treated perturbatively. 
In particular, the pairing interaction can be chosen to 
be a function of the quasiparticlc spins and momenta 
near the Fermi surface unaffected by the SO coupling. 
This approach to the theory of noncentrosymmetric su- 
perconductivity was introduced by Edelstein in Ref. 7 
and further developed in Refs. 8.9i lTolllTI . Due to the 
absence of inversion symmetry the superconducting or- 
der parameter does not, in general, have definite parity, 
becoming a mixture of spin-singlet and spin-triplet com- 
ponents. The triplet component appears, however, only 
when a spin-triplet channel is explicitly present in the 
pairing interaction, even for a finite value of the SO band 
splitting. When the SO coupling increases, only a certain 
type of the triplet pairing survives, along with the singlet 
component. 

In the limit of strong SO coupling, i.e. when the band 
splitting exceeds the superconducting energy scales, the 
Cooper pairing between the electrons with opposite mo- 
menta occurs only if they are from the same nondegen- 
erate band. Interband pairing is still possible, but only 
at the band degeneracy lines or points in the momen- 
tum space, and can be neglected. According to Ref. [H, 
this is the scenario that is realized in CePtaSi: The band 
structure calculations show that Ego ranges from 500K 



to 2000K, which is much larger than the critical tem- 
perature T c = 0.75K, even taking into account a possi- 
ble mass renormalization due to the strong correlation 
effects. The same is likely to be the case in other mate- 
rials, e.g. L^PdsB and L^PtaB, see Ref. [l3|. The pair- 
ing interaction in the strong SO coupling case is most 
naturally introduced using the basis of the exact band 
states , 12 i 14 i 15 i 16 which already incorporate all the effects 
of the crystal lattice potential and the SO coupling. 

In the band representation, the superconducting or- 
der parameter is given by a set of complex gap functions, 
one for each band, which arc coupled due to the interband 
scattering of the Cooper pairs and other mechanisms, e.g. 
impurity scattering. The overall structure of the the- 
ory resembles that of multi-band superconductors ! 17 > 18 
However, since the bands are nondegenerate, the pair- 
ing symmetry is peculiar: While each order parameter 
is an odd function of momentum, the gap symmetry, 
in particular the positions of the nodes, is determined 
by one of the even representations of the point group of 
the crystal. When expressed in the spin representation, 
the order parameter becomes a mixture of singlet and 
triplet components^ the latter appearing even without 
any spin-triplet term in the pairing interaction, as an in- 
evitable consequence of the SO band splitting and the 
difference between the gap magnitudes and the densities 
of states in different bands. 

Noncentrosymmetric superconductors exhibit a variety 
of unusual features, which are absent in the centrosym- 
metric case, such as a strongly anisotropic spin suscepti- 
bility with a large residual component,— d°i 14 i 19 i 20 i 21 i 22 i 23 
magnetoelectric effect,— & & 20 i 24 i 25 i 26 an d helical super- 
conducting states 2Z£&22£±£Li££2. Other properties of in- 
terest include the nuclear spin relaxation , 34 ! 35 Joseph- 
son and quasiparticlc tunnelling, 36 - 37 electron correlation 
effects , 25 ! 38 superfluid density and the London penetra- 
tion depthj 3 ^ and the free energy in the clean case^i 

In this article we study the effects of nonmagnetic im- 
purities on the superconducting properties. In contrast 
to the treatment of the same problem in Refs. [ll| and 
l2tl we assume that the SO coupling is larger than the 
superconducting energy scales, which necessitates using 
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the band representation of the pairing Hamiltonian. Our 
main goal is to find how the superconducting critical tem- 
perature and the upper critical field depend on the im- 
purity concentration, for both conventional and uncon- 
ventional pairing. The article is organized as follows. In 
Sec. HI we obtain the band representations of the im- 
purity Hamiltonian and of the disorder-averaged Green's 
function of electrons in the normal state. In Sec. IIII1 we 
derive the Ginzburg-Landau (GL) free energy functional 
in the superconducting state. In Sec. IIV1 we calculate 
the critical temperature and the upper critical field in the 
presence of impurities and also discuss the application of 
our results to a model of superconductivity in CePtsSi. 



II. IMPURITY SCATTERING IN NORMAL 
STATE 

In a noncentrosymmetric crystal with SO coupling the 
electron bands are nondegenerate. The formal reason 
is that without the inversion operation one cannot, in 
general, have two orthogonal degenerate Bloch states 
at the same wave vector k. In the limit of zero SO 
coupling there is an additional symmetry in the system 
- the invariance with respect to arbitrary rotations in 
spin space - which preserves two-fold degeneracy of the 
bands. Let us consider one spin-degenerate band with 
the dispersion given by eo(fe), and turn on the SO cou- 
pling. The Hamiltonian of noninteracting electrons in the 
presence of scalar impurities can be written in the form 
H = H a + H imp , where 

H o = ^2 [eo(k)5 a p + -/(k)cr al 3}al a a k f3, (1) 

fc,a/3 

a, (3 =T, I denote the spin projection on the z-axis, J^k 
stands for the summation over the first Brillouin zone, a 
are the Pauli matrices, the chemical potential is included 
in eo(fc), and 

H imp = f d 3 r£tf(r)V4(rhMr). (2) 

The impurity potential U(r) is a random function 
with zero mean and the correlator (U{ri)U(r2)) = 
nimpU§5(ri — T2), where rii mp is the impurity concen- 
tration, and Uq is the strength of an individual point-like 
impurity. 

The "bare" band dispersion satisfies eo(—k) = e (k), 
eo(g~ 1 k) = eo(k), where g is any operation from the 
point group G of the crystal. The SO coupling of elec- 
trons with the crystal lattice is described by the pseu- 
dovector 'y(fc), which has the following symmetry proper- 
ties: -f{k) = — *y(— k), (g-f){g~ 1 k) = j(k). For example, 
the point symmetry of CePtaSi, CeRhSi3 and CeIrSi3 is 
tetragonal and described by G = C4,,, which is generated 
by the rotations C^ z about the z axis by an angle 7r/2 
and the reflections a x in the vertical plane (100). The 



pseudovector -f(k) can be written as 

7 (fc) = 7_l [<j)E,u (fc) X Z] +7\\<t>A 2 ,u(k)z, (3) 

where j± and 711 are constants, and 4>e,u and 4>a 2 ,u are 
the odd basis functions of the irreducible representations 
E (two-dimensional) and Ai (one-dimensional) of C4,, ^ 
The Hamiltonian ([1]) with ~f(k) given by Eq. {3]) is a 
three-dimensional generalization of the Rashba model, 
which is widely used to describe the effects of SO coupling 
in two-dimensional semiconductor heterostructures4^ 

The SO coupling strength depends on the quasipar- 
ticle momentum and might vanish, for symmetry rea- 
sons, along some directions or at some isolated points 
in the Brillouin zone. The former possibility is real- 
ized in the tetragonal case: 7(fc) = along the axis 
k x = k y = 0, which can be seen from the polynomial ex- 
pressions for the basis functions: <pE,u(k) ~ (k x ,k y ) and 
4'A2,u(k) ~ k x k y k z (k^, — ky). However, this is not generic: 
For example, the point symmetry of Li2(Pdi_ a; ,Pt :I ;)3B is 
described by the cubic group G = O, which contains only 
the rotations about the axes of the second, third, and 
fourth order. In this case the function f(k) transforms 
according to the vector representation F± : 

7(fc) =7o0F 1 ,«(fc). (4) 

The representative expression for the basis function is 
simply <j>p 1<u (k) ~ (k x ,k y ,k z ), so that the SO coupling 
vanishes at the point k x = k y = k z = 0. 

The diagonalization of Hq yields the following eigen- 
values and eigenstates: 

£ A (fc) = eo(fc)+A| 7 (fc)|, (5) 

and 

Xkx{r,a) = ^=u aX {k)e lkr , (6) 
where A = ± is the band index, V is the system volume, 

v/2| 7 |(l7l + A7*) 

and 9\ are arbitrary (in general fc-dependent) phases. 
The Bloch spinor components u a \ form a unitary ma- 
trix u(k). It follows from Eq. ([5]) that in the presence 
of SO coupling the degeneracy of the electron bands is 
lifted everywhere in the Brillouin zone, except maybe for 
some high-symmetry lines or points, where j(k) — (the 
band structure might be such that the zeros of 7 are 
not located on the Fermi surface). For the typical band 
splitting Ego one can use j f° r instance, the Fermi-surface 
average of 2|— | . 

The band representation of the free-electron Hamilto- 
nian has the following form: 

ffo = EE uk)ci x c kX . (8) 

k A=± 
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The band dispersion functions ([5]) are even in k due to 
time reversal symmetry: the states \kx and K\kx belong 
to k and —k, respectively, and have the same energy. 
Here K = i&2Ko is the time reversal operation, and Kq 
is the complex conjugation. 

Writing the field operators in Eq. |(2J) in the form 
''Pair) = J2k,\Xk\ir,a)c k \, we obtain the band repre- 
sentation of the impurity Hamiltonian: 



kX Ck>\> 



(9) 



"4- 



x 3 x 4 



kk' XX' 



FIG. 1: The impurity self-energy in the band representation. 
The dashed line corresponds to ntmpUQ, the vertices include 
the anisotropy factors w(k,k'), and the solid line is the av- 
erage Green's function of electrons in the normal state. It 
is shown in the text that the self-energy is nonzero only if 
Ai = A4 and A2 = A3. 



where 



Uxx'ik,k') = Uik-k')wxx'ik,k'), 



(10) 



U(q) is the Fourier transform of the impurity potential, 
{U(q)U(q')) =n lmp U$V6 q - ql , and 



w(k,k') = u\k)u(k'). 



(11) 



We see that the impurity scattering amplitude acquires 
both intraband and interband contributions and also be- 
comes anisotropic, even for isotropic scalar impurities. In 
more general cases, for instance in the presence of spin- 
dependent scattering in Eq. @ , the band representation 
of the impurity Hamiltonian still has the form ([9]), but 
with a different Uxx'ik,k'). 

We note that the matrix elements of the anisotropy 
factor w satisfy the following useful identity: 



{wxx'i^k'f 



1 + XX'j(k)j(k') 



(12) 



which can be easily verified using the explicit expressions 
([7]) for the Bloch spinors. 

The electron Green's function in the band representa- 
tion is introduced in the standard fashion^ 



Gxx'ik,T;k',T') = -<T T c feA (r)4 A ,(T')>. 



(13) 



In the absence of impurities it has the following form 

Sxx' 



G 



O.AA 



UkY 



(14) 



where uj n = (2n + l)7rT is the fermionic Matsubara fre- 
quency (we use the units in which ks = 1). Now we 
will show that, even when the disorder is included, the 
average Green's function remains band-diagonal. 

The disorder averaging with the Hamiltonian ^ can 
be performed using the standard methods, see Ref. Hl|. 
The result is the matrix Dyson equation of the form 



G- 



G n 1 



E, where G is the average Green's function 



and £ is the impurity self-energy. In the Born approxi- 
mation, 



£(fc,cj„) = n 
d 3 k' 



U 2 

imp u 



(27T) 



w(k,k')G(k',uj n )w(k',k) (15) 



(here we have taken the thermodynamic limit V — * 00). 
The diagrammatic representation of the self-energy is 
given in Fig. [TJ We seek solution of the Dyson equa- 
tion in a band-diagonal form: G\\> — G\6\\>, then the 
integrand on the right-hand side of Eq. (|15p becomes 

u{k')G{k' \u n )ii){k') 

G + (fe',^ n ) + G_(fe', Wn )„ 
= 2 To 

+ G+(fc> " ) ^- (fc> " ) C(fcQ i (16) 

where Ti are the Pauli matrices, and £(fc) = u(k~)T^u^ '(k) . 
Using the expressions (|7|), one obtains £(fc) = 7(fc)r, 
which is an odd function of k (7 = 7/I7I). Therefore 
the last line in Eq. (|16p vanishes after the momentum 
integration, and T,(k,uj n ) = E(w n )fo. The real part of 
the self-energy renormalizes the chemical potential, and 
for the imaginary part we obtain: 



Im£(cj n ) = —Tmim P U 



N+ + N- 



■ signup 



where N x = V -1 ^2 k 5[£\(k)] is the Fermi-level density 
of states in the Ath band. In order to preserve the usual 
form of the impurity self-energy in the normal state, we 
introduce the notation Np = (N + + N-)/2 and also de- 
fine the elastic mean free time as r = {2^ni mv Ul Nf)~ x ■ 
Then the average Green's function takes the form 



G\x'{k,uj„) 



'XX' 



iu n - £,x(k) + (i/2r) signw r , 



(17) 



The derivation above is valid under the assumption 
that the elastic scattering rate is small compared with the 
Fermi energy , which justifies neglecting the diagrams 
with crossed impurity lines in the self-energy in Fig. [T] 



A. Effects of magnetic field 

In the presence of a nonzero uniform magnetic field B 
Eq. is replaced by 



H = ^2ci x Sxik) 



Ckx, 



(18) 



k.X 
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where £\ is the effective band Hamiltonian in the k- 
space42 It was first pointed out by Peierls^ based on the 
requirement of gauge invariance, that the magnetic field 
can be included in the band electron theory by simply 
replacing the wave vector k in the zero-field band disper- 
sion by the operator K = k + (e / he) A(f) , where 
r = iVfc is the position operator in the fc-representation 
and e is the absolute value of the electron charge. We 
use the symmetric gauge, in which 

K = k + i^-(BxV k ). 

2hc 

The Peierls substitution gives only the zero-order term 
in the expansion of the effective single-band Hamiltonian 
in powers of B, which in our case can be written as 



£ x (k) = UK) - Bm x (K) 



(19) 



The second term on the right-hand side is the analog of 
the Zeeman interaction for nondegenerate bands. ?2 In the 
generalized Rashba model ([TJ) it has the following form: 



m\(k) = XfiBj(k) 



(20) 



(/its is the Bohr magneton), which is valid everywhere 
except for the vicinity of the band crossing points, where 
the approximation of independent nondegenerate bands 
fails and the effective band Hamiltonian approach might 
not work. 

It is convenient to introduce the Fourier transform of 
the Green's function CC; 



Gxx>{r,r';uj n ) = I ^ e lfcMfcV G A v(fc, k';u n ). (21) 



kk' 



We would like to stress that this is not the same as 
the electron Green's function in the coordinate-spin 
representation. The latter is defined as (ra\(iu) n — 
H)- l \r'(3) = Z kk , tXX ,{ra\k\)G xx ,(k,k';oj n )(k'\'\r>P), 
where (ra|feA) = Xk\(v,a) is the spinor wave function 

©, witha,/3=t4- 

The Green's function (f2"Tj) satisfies the equation 

(iw n r -ho - him P )G(r, r'\ u„) = f Q S(r - r'). (22) 

Here ho,x\> — S X x'£x(K r ) is obtained by replacing K in 
Eq. Gjjj by 



K r = -iV r + -^-A(r)=-iV r + 7 ^-(Bxr), (23) 
he 2hc 

and hi m p is the Fourier transform of the impurity Hamil- 
tonian ([9]). The latter is nonlocal in real space and can 
be represented as a differential operator of infinite order: 

{him P G)xx'{r,r';uj n ) =^[W 0iAAl (r) 

Ai 

+U liXXl (r)K r + ...]G Xl x(r,r';u n ), (24) 



where 

U h xx>(r) 



1 E' 



v 



^(f.-f), 



V 



- V e ipr 

V £—1 

P 

x2(V p -^ 



U X \> 



p p_ 

2' 2 



etc. The magnetic field affects the electron Green's func- 
tion in Eq. |22|) through the vector potential in the op- 
erators K r (the "orbital effect") as well as through the 
m-term in Eq. (fl9|) (the "paramagnetic effect" ) . 

In the absence of impurities the solution of Eq. p2|) 
is band-diagonal and can be represented in a factorized 
formji^ 



G X X'(r,r']uj n ) = G xx > (r , r' , w n ) 

e 



x exp 



i— I A(r)dr 

he 



, (25) 



where the integration in the phase factor is per- 
formed along a straight line connecting r and r', and 
G xx >(r,r',oj n ) = 5xx'G ,x(r - r',u> n ). In_superconduc- 
tors the orbital effect of magnetic field on Go can usually 
be neglected^ and one obtains: 



G .x(k, lu„ 



Uk)+m x (k)B 7 



(26) 



see also Ref. |3l|. 

The same argument can be used also in the disor- 
dered case, to show that the Green's function before 
disorder averaging has the form ([2^)1 . where Gxx' sat- 
isfies Eq. (f2"2"|) . in which the vector potential is formally 
set to zero, but the paramagnetic term is still present. 
In order to perform disorder averaging of Gxx 1 , we go 
back into the band-momentum representation (|13[) and 
repeat the steps leading to Eq. (fT5")) . The only differ- 
ence is that now the Green's function is no longer even 
in k due to the presence of the paramagnetic term with 
m x (k) = —mx(—k), see Eq. (|2"6")) . Therefore the last 
line in Eq. (|16p does not vanish identically after the 
momentum integration, and the impurity self-energy is 
no longer band-diagonal. One can show however that 
the off-diagonal contributions to £ A a' are °f the order of 
fisBN'p, which can be neglected. Therefore we arrive at 
the following expression for the impurity-averaged Gxx' '■ 



(Gxx> (k, k'; Lo n )) imp = Sxx>5 k ,k'Gx{k, u n ), (27) 



where 



G x (k,uj n ) 



■£\(k) +m x (k)B + (i/2r) signw„ 



(28) 



While the orbital effect of the field is reduced to the phase 
factor in the Green's function, see Eq. (|2l>)) . the Zeeman 
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interaction survives in G\y and affects the quasiparticle 
energies. The Zeeman term will be eventually dropped, 
because its effect on the GL free energy turns out to be 
negligible, see the next section. 



III. SUPERCONDUCTING STATE 

Now let us take into account an attractive interaction 
between electrons in the Cooper channel, using the basis 
of the exact eigenstates of the noninteracting problem. 
The total Hamiltonian is given by H = H +H imp + H int , 
where the first two terms are given by Eqs. © and 
and the last term has the following form: 

Hint = 2V ^ ^ Vxx '( k ' fe/ ) c l+<j/2,A c -fc+q/2,A 

kk'q XX' 

'XC_ k < +q /2,X'Ck'+q/2,X>- (29) 

We assume, in the spirit of the Bardeen-Cooper-Schrieffcr 
(BCS) theory, that the pairing interaction is nonzero only 
inside the thin shells of width e c in the vicinity of the 
Fermi surfaces, i.e. when |^(fe)|, |£ A /(fc')| < e c . The 
pairing potential satisfies V\\i(— k,k') = V\y{k,—k') = 
—V\x>(k,k'), which follows from the anti-commutation 
of the fermionic operators. The diagonal elements of the 
matrix V describe the intraband Cooper pairing, while 
the off-diagonal ones correspond to the pair scattering 
from one band to the other. The SO splitting of the bands 
is assumed to be large compared with all the energy scales 
associated with superconductivity. In this case the for- 
mation of the pairs of electrons belonging to different 
bands is strongly suppressed. Although the bands may 
touch at some isolated points at the Fermi surface, the 
interband pairing in the vicinity of those points is still 
suppressed due to the phase space limitations. 

The pairing potential can be represented in a factorized 
form: 

Vx\>{k,k') = -Vxx>tx(k)t* y (k')MkWx>{k'), (30) 

where the coupling constants Vxx' form a symmetric 
positive-definite 2x2 matrix, t\(k) = —t\(—k) are 
non-trivial phase factors, see Refs. Il4ltl5l and <^a(&) = 
^Tg (k) are even basis functions of an irreducible repre- 
sentation r of the point group G of the crystal (we keep 
only the representation which corresponds to the pair- 
ing channel with the maximum critical temperature) 4£ 
While 4>+(k) and (f>-(k) have the same symmetry, their 
momentum dependence docs not have to be exactly the 
same. We consider only one-dimensional representations 
of G, because, to the best of our knowledge, there is 
no experimental evidence of multi-dimensional supercon- 
ductivity in real noncentrosymmetric materials. 

The basis functions are nonzero only inside the BCS 
shells and are normalized: (\cj)\(k)\ 2 )\ = 1, where the an- 
gular brackets denote the averaging over the Fermi sur- 



face in the Ath band: 

(("OK = —£(■..)«*=)]• 
k 

The phase factors t\(k) originate in the expression for the 
time reversal operation for nondegenerate band electrons: 

K\k\) = t x (k)\ -k,X). (31) 
For instance, for the eigenstates (O we obtain: 

t x (k) = ^-*[M»)+M-»)1 7 x(fc)-Hy(fc) , (32) 

Since the factors t\ explicitly depend on the arbitrary 
phases 9\(k), they must drop out of the final expressions 
for all observable quantities. 

It is instructive to see how the phase factors t\(k) ap- 
pear in a simple BCS-like model in which the pairing 
interaction is local in real space: 

#*»t = -T E [ d 3 r^ a (r)^l(r)Mr)Mr), (33) 

"/3=U 

where V > is the coupling constant. Using the band 
representation of the field operators, we obtain: 

ff,n ' = 2V E E t/ A 1 A 2 A 3 A 4 (fc,fc') 
kk'q Ai, 2,3,4 

XC fe+q/2,A 1 C -fe+q/2,A 2 C -fe'+'?/2.A 3 C fc'+q/2,A4 7 (34) 

where 

UA 1 A 2 A 3 A 4 (fe,fc') = -V^«* Al (fe)^ Aa (-fe) 

a/3 

xuf3x 3 (-k')u a \ 4 (&'). 

Here we replaced u a \(±k + q/2) by u a \(±k), neglecting 
the corrections of the order of 0(q/kp). The expression 
(j34[) contains both intra- and interband pairing terms. 
For the reasons explained above, we neglect the latter 
and set Ai = A2 = A and A3 = A4 = A', which reduces 
the Hamiltonian to the form (|29|) . Using the identities 

u a x{-k) = t* x {k) J2( icr 2) a /}U* 0X (k), (35) 

P 

which follow from the definition (|3"Tj) of the phase factor 
t A (fe), we obtain: 

Vxy(k,k') = -Vtx(k)t* x ,(k')\wxy(k,k')\ 2 

see Eq. p^|) . The terms in V containing the 7's are even 
in both k and k' . Using the anti-commutation of the 
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fermionic operators one can easily show that these terms 
drop out of the Hamiltonian. Finally, 

tf« = -7y£E*A(*0*A'(fc') 

kk'q AA' 

XC l+q/2,A C -fc+q/2,A C -fc'+q/2,A'Cfc'+q/2,A'-(36) 

Thus the band representation of the pairing potential 
in the model (|3"3")) necessarily contains the phase factors 
t\(k). The gap symmetry corresponds to the unity rep- 
resentation with 4>\{k) = 1, and all coupling constants 
take the same value: V\\> = V/2. 



k+q/2 k+q/2 

= \ 

-k+q/2 -k+q/2 



k'+q/2 




-k+q/2 -k'+q/2 



FIG. 2: The diagrammatic representation of the gap equation. 
The triangle denotes the gap function Ax, the filled triangle 
- the impurity-renormalized gap function D\, the dashed line 
corresponds to W\\i (k, k'), and the solid lines are the average 
Green's functions of electrons. 



A. Derivation of the free energy functional 

The superconducting order parameter in the static case 
can be represented as 



A x (k,q) =r}x(q)t x (k)(f>x(k), 



(37) 



where the bosonic fields rj\ play the role of the order pa- 
rameter components. The free energy T (more precisely, 
the difference between the free energies of the supercon- 
ducting and the normal states at the same temperature 
and field) is a functional of rj\. In the vicinity of the tran- 
sition temperature at arbitrary field the order parameter 
is small, and one can keep in the free energy expansion 
only the terms quadratic in rjx- Following the proce- 
dure outlined, e.g. in Ref. l3ll . we obtain the impurity- 
averaged free energy in the form T = T\ + Ti, where 



* = JE^E"A(«)^v»?v(g), 



2 ^ V 

AA' 



(38) 



and 



* = -^E T E^E A A(^)Av(fc', 9 ) 

AA' n kk'q 

x<5 AA ,(-fc+|,-fc'+|;-u; r 



imp 



(39) 



This expression is applicable for any pairing symmetry 
described by a one-dimensional representation of G, and 
for any band structure. There are corrections to T% of 
the order of (T c /e F ) 2 , related to the orbital magnetism of 
the Cooper pairS ) 31 i 46 which we neglect. The GL gradient 
expansion in all orders can be obtained from the Taylor 
expansion of Eqs. (|3"8|) and (f3"9"| in powers of q, by making 
the replacement 



D = -i\7 r 



2e 

the 



A(r) 



(40) 



in the final expressions. 

The disorder averaging in Eq. (|39[) involves sum- 
ming the impurity ladder diagrams, in which the aver- 
age Green's function is given by Eqs. ([2T]) . (|2"8)l . In the 



thermodynamic limit V — > 00, the momentum sums are 
replaced by integrals, and we obtain: 



d 3 



d 3 k 



ly . , 7 , 

(27r )3 Z^J (27r) 3 

xD x {k,q,uj n ), (41) 



r A A (fe,g) 



-k + 2' _CJ « 



where Dx are the impurity-renormalized gap functions, 
which satisfy the following integral equations: 



D x (k,q,uj n ) = A x {k,q) + y 



j^Wxx.{k,k>) 



(V v (// + |,u, n )G A ,(-fc' + |,- Wn 



with 

Wxx'(k,k') 
1 



1 



xD x ,(k',q,u n ),(42) 
wxx> (k, k')wxx> {-k, -k') 



2ttN f t 



2ttN f t 

tx(k)t$,(k')wxx'(k,k')wx'x(k',k) 



1 



2ttN f t 



h(m ^ } i±^mi (43) 



corresponding to the impurity line, see Fig. [2J To obtain 
the last expression we neglected the difference between 
w(k + q/2, k' + q/2) and w(k, k') (this is legitimate be- 
cause the order parameter varies on the scale q <C kp , 
and w is a smooth function of fc, k' near the Fermi sur- 
face), and also used the identities (j35]) and the fact that 
w^(k',k) = w(k,k'). 

We seek solution of Eq. |42|) in the form 

Dx(k,q,co n ) = tx(k)Ax(k,q,co n ). (44) 

Due to the momentum cutoff built in the basis functions, 
A A is restricted to the vicinity of the Ath Fermi surface. 
Since the gap functions weakly depend on energy in the 
vicinity of the Fermi surface, one can integrate the prod- 
ucts of two Green's functions with respect to £ A = £,x{k): 



k + 



Nx I d£xGx(k + i,u; n )Gx 



irN F Lx(k,q,u n ), 
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where 

L x (k,q,u) n ) 



ftx 



\u) n \ + 1/2t + iCl X (k, q) signw„ 



(45) 



depends on the direction of k. The notations are as fol- 
lows: 



P\ 



(46) 



is the fractional density of electronic states in the Ath 
band (0 < p A < 2, p+ + p- = 2), 

n x (k,q) = ^P^-m x (k)B, (47) 

and v\(k) — hr 1 d^ x (k)/dk is the Fermi velocity. In this 
way we obtain the following equation for A a as a function 
of the direction of k, at given q and lu u : 



A x (k,q,u n ) = rix{q)(j)\{k) 
1 

+ 2~r 

xA X '(k',q,uj n ) 



a: M 9 Ly(k,q,uj n ) 



(48) 



The phase factors £a(&) drop out of the expression (j4Tj) 
for J^2, which takes the following form: 



nN F f d 3 



xTW^(fc)L A (fc ) «i J w n )A A (fc,q,« n )) . (49) 



We note that only the diagrams containing G+G+ and 
G_G_ have been included in the impurity ladder in Eq. 
([42]) . The contribution from the diagrams with G+G_ 
is smaller by a factor of max(e c , t~ 1 )/Eso "C 1 and ne- 
glected. 

The equations for A x (k, q, w n ) can be solved by trans- 
forming them into a system of linear algebraic equations. 
We introduce 

X a (q,oj n ) 

= (A x , a (k)L x (k, q, w n )A A (fe, q, lo 71 )) x , (50) 

A 

where a = 0, 1, 2, 3, and 

A (h\ _ / 1 , a = 

AA,aW-< A%(fe) ? a = 1)2 ,3 ■ 



Then Eqs. (|48p can be written as 

3 

4r 



1 3 

A A = ?7A0A + -7- ^2 ^X a,X a , 



(51) 



and the contribution (|49f to the free energy takes the 
following form: 



T 2 



N F f d 3 q 



2 J (2tt) 3 

^b ? A( 9 )|Vr^(|0A(fc)| 2 iA(fe,g,c„)) ; 

A n 



A n a— 

x {^ A (fe)A AiQ (fe)i A (fe,g,w n )) ; 



(52) 



From Eqs. (|50|) and (T5"l"j) we obtain a system of four 
algebraic equations for X a (q,u) n ): 



J2{S ab - B ab )X b = Y a , 



(53) 



6=0 



where 



B a b(q, w n ) = j- ^2 (A x . a (k)A x ^ b (k)L x (k, q, v n )) x , 

A 

Y a (q, uj n ) = Vx(q) (<}) X (k)A Ka (k)L x (k, q, w n )) A . 



We see that, since Y a depend linearly on 77A, so do 
X a , and the right-hand side of Eq. ([52")) is a bilinear 
functional of the order parameter components. The total 
free energy, which also includes the contribution (|38p , can 
be written in the following form: 



T = \ ?|^3-X]^(g)/AA'(g,-B)f7A'(g) 



(54) 



The explicit expressions for the coefficients /aa' can be 
derived in principle at arbitrary q and B, using the solu- 
tion of the linear equations (|53|) . This procedure would 
allow one to calculate the upper critical field H C 2 at any 
temperature. We focus on the weak-field limit in the 
vicinity of the critical temperature, in which case Eqs. 
([53f can be solved using a gradient expansion, see Sec. 
HUB] below. 

We note that the free energy of our system is formally 
equivalent to that of a two-band singlet superconductor, 
in which the impurity scattering amplitude, the pairing 
symmetry, and the Zeeman coupling are all anisotropic. 
Significant simplifications are achieved, for example, in 
the case when the SO band splitting is so large that 
there is only one nondegenerate band (say, the "+"-band) 
crossing the Fermi level. This is formally described by 

setting p_ = and V-\ = V-- = 0. We will refer to 

this case as the single-band limit. On the other hand, in 
the model (|33[) the matrix V is degenerate and the whole 
formalism should be modified. One can show that in this 
case the pairing interaction l|36p can be decoupled using 
just one bosonic field 77(g), so that the gap functions are 
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given by A\(k,q) = i](q)t\(k), and the GL theory has 
the same form as in the standard, i.e. one-component 
isotropic singlet, case. The physical properties of these 
special cases will be studied below by taking the appro- 
priate limits in the general two-band expressions. 



B. Gradient expansion 

The gradient and field expansion of the free energy (154)) 
is obtained by expanding both L\ and X a in Eq. (|52| in 
powers of Q\. In order to solve Eq. (|53p it is convenient 
to represent it in the symbolic form (1 — B)X = Y. The 
4x4 matrix B can be written as B — Bq+SB, where in the 
first term £l\ is set to zero, while the second term contains 
the Independent corrections. Similarly, Y = Y n + SY 
and X = X Q + SX, where X = (1 - Bq^Yq. Then, 

SX = (1-Bq -5B)-\5Y + SBX ) 
= {l-B^^Y + SBX ) 

+ (1 - B )- 1 5B{1 - Bo^iSY + 5BX ) + ... 

Inserting this in Eq. (jB"2")) , we obtain after some lengthy 
but straightforward algebra: 



/ = fa + h + 



(55) 



where f m denotes the terms of the order of f2™. The 
terms with odd powers of £l\ are proportional to sign u> n 
and vanish after the summation over the Matsubara fre- 
quencies. 

The uniform and field-independent contribution is 
given by 



A\\> — fo.xX' — -r— V xx ) — p\5\\>S l 

l LIXp 



'-rP\py{4>\){4>\')Sn 

4t 



(56) 



To make the notations compact, here and below we omit 
the arguments of the basis functions and the subscripts A 
in the Fermi-surface averages, and also assume that the 
basis functions are real. The Matsubara sums $k,i are 
defined as follows: 



i 



LO n \ k {\LO n \+l/2T) 1 ' 



(57) 



The logarithmically divergent sum in Sq\ is cut off at 

n c ~ e c /2nT > 1: 



Soi =2*TJ2 



1 



n=0 



-1/2t ~ ttT 
Here In 7 ~ 0.577 is Euler's constant, 

F(x) 



= ln^-F(rT). (58) 



V 2 Aizx ) 



(59) 



and 9(x) is the digamma function?^. One can also check 
that S u = 2t¥(tT). Therefore, we have 



.4 



XX' 



N F 



1 T/-1 A 1 2 ^ £c 

—V xx ,-px6xx'ln — 



+ {pxSxx- - -paPA'<0a)<^))f(tT)]. (60) 

The second-order term in the expansion (|55[) has the 
following form: 



where 

Txx>(q,B) 



: ((0 A )(^r! 2 A ,) + (0v)(0A^»^l 

Ai 

-AA' Y, (%^x)(^x'^x')Si 



(61) 



03./ 



i— x,y,z 



+i E (Ew*n*>) 

i—x,y,z Ai 

X {\'{<px){li<t>X>nX>) + A(0A')(7*0Ar>A» §13,i 
1 ,,,,,, . ~ ,\ 2 ^ 



(4r) 2 



i=x,y,z Ai 



with 



Ski.i = 7rT V] - — iju 71 — n — r 
^ |w„| fc k + 1 



\to n \ k (|w„| + 1/2t)' |w„| + Si /2t 



, (62) 



and Si = l--EAPAa 2 (fe))A/2. 

Along with the usual second-order gradient terms, the 
free energy expansion also contains the terms linear in 
both q and B, which give rise to the magnetoelectric ef- 
fect, see Refs. SHU], and also lead to the helical super- 
conducting phases considered in Refs. 29l3lll32l It can 
be shown, however, that these effects are small. Indeed, 
a typical contribution to f% is of the form 7V^(i7 A )/T 2 . 
Minimizing it with respect to q, we find q ~ pbH/Hvf, 
therefore the correction to the energy due to the Zee- 
man term in Eq. (|47| is proportional to (pbH/T c ) 2 ~ 
(T c /eF) 2 (m* /m) 2 (H/H C 2) 2 , where m* is the effective 
mass of quasiparticles and H C 2 is the upper critical field 
at T = 0. This ratio is typically very small (unless the 
smallness of T c /ep is compensated by a large value of 
the effective mass), so it is legitimate to drop the Zee- 
man terms and use 



Qx(k,g) 



hv x {k)q 



(63) 



As a result, f 2 becomes a second-degree polynomial in 
q- f2,XX'(q) = J2ij K\\' t ijqiqj, where the coefficients arc 
found from Eq. (|6ip . 
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Finally, returning to real space and making the sub- 
stitution q — > D, see Eq. (|40p . we obtain the GL free 
energy in the second order of the gradient expansion: 



A 



XX' 



J2 K 



xx',ijDiDj 



Vx'(r), 



(64) 

with i,j = x, y, z. The explicit expressions for the expan- 
sion coefficients are determined by the pairing symmetry 
and given below. 

Unconventional pairing: The order parameter trans- 
forms according to a nontrivial one-dimensional represen- 
tation of the point group of the system, and ((j)\(k))\ — 
0. From Eq. (pO)) we obtain: 



Axx> = 



1 t/-i A i 2 7 e c 
W V xw ~ P x5xx' In 



AA , ... ^ 

+ P xSxy¥{ T T) 

while many of the terms in Eq. (|61[) vanish, giving 
h 2 N F 



(65) 



K 



XX', ij 

1 
17 



\pxSxx'{<Px v x,i v x,j)So3 

k<t>x>vx',j)S 3,k ■ (66) 



k—x,y,z 



Conventional pairing: The order parameter transforms 
according to the trivial (unity) representation of the 
point group. For simplicity, we consider only a fully 
isotropic pairing, for which <fi\(k) — 1, then 



Axx> = 



j^V xx ) - PX 5 XX , In 



2t£c 
ttT 



and 



Kxx',ij = 



h 2 N F 



~ P+P _¥(tT) 



pX§XX'(vx,iVX,j)So3 



(67) 



1 



{(vx,iVX,j) + (vx',iVX',j})Si3 



+ ^-(P+( v +,*v+. 3 ) + p-(v- ti v- tj ))S 2 3 
+AA' ^ (%vx,i)(%vx',j)So3,k 

k—x,y,z 

+ i7 (P+(^kV+,i) - P-(%V-,i)) 

k=x,y,z 

x(M%Vx,j) + X'(%Vx',j))Si3,k 

1 



(4r)2 



(68) 



fc— x,y,z 

x(p + (%v +J ) - p-(j k v-j))S 2 3,k 



In the next section we use the GL expansion (|64|) to 
study the effects of impurities on the critical temperature 
and the upper critical field. 



IV. CALCULATION OF OBSERVABLES 
A. Critical temperature 

At high temperatures the matrix A in Eq. (|64[) is pos- 
itive definite, therefore the minimum of the free energy 
is achieved at 77+ = r\- = 0, which corresponds to the 
normal state. The superconducting critical temperature 
T c , defined as the temperature below which one of eigen- 
values of A turns negative at zero field, is found from the 
equation det A = 0. 

It is convenient to introduce the following notation: 



gx x' = N F V\x>p x > = Vxx'Nx> 



(69) 



Since V is positive definite, we have g++ > 0, (/__ > 0, 
det <7 > (note that the matrix g is not symmetric, in 
general). Using the expression (jBT))) . the equation for the 
critical temperature can be written as det(fo + gM) = 0, 
where 



M 



XX' 



JXX' 



hi 



2t£c 

7tT c 

Px 



+ {Sxx'-^f{4>x){4>x'))nrT c ). (70) 

At r = 00 the second term in M X x> vanishes, and we 
obtain the critical temperature of the clean superconduc- 
tor: 



where 




(71) 



(72) 



is the effective coupling constant. 

We note that in the model (|33p , in which the pairing is 
described by a single coupling constant Vxx' — V/2, we 
have g = NpV. Although the expression for the critical 
temperature in this case has the usual BCS form, the su- 
perconductivity is non-BCS, because the order parameter 
resides in two nondegenerate bands (see the discussion in 
the end of Sec. IIII A[) . with the critical temperature in- 
dependent of the band splitting. 

In the presence of impurities the cases of conventional 
and unconventional pairing have to be considered sepa- 
rately. 

Unconventional pairing: = 0. In this case, we 

obtain the following equation for T c : 



ln^=F( T T c ), 



(73) 



where F(x) is defined by the expression (|59|) . The reduc- 
tion of the critical temperature is described by a universal 
function, similar to the suppression of sup erconductivity 
by paramagnetic impurities, see Ref. 48. In particular, 
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at weak disorder, i.e. in the limit tT c q 3> 1, we have 
F(tT c ) ~ it/8tT c0 , therefore 



T c — T c q — 



Using the small- x asymptotics 



7T 

8r' 



(74) 



W(x) = In ( -L) + ^x 2 + 0(x 3 ), (75) 



we find from Eq. (|T3|) that the superconductivity is com- 
pletely suppressed at tT c q = j/ir ~ 0.567. 

Conventional pairing: cf>\ — 1. Instead of Eq. (]T3[) we 
obtain: 



In- 



r, 



l + Cl F(x) 



1 



c 2 + c 3 F(x) + v /c 4 + c 5 F(x) +c 6 F 2 (x) 5 
where x = tT c , and 



--, (76) 



Cl = 


P+(.9 


— - 9+- 






C2 = 


.9++ 


+ 9— 




2 


C3 = 


dot g 
2 




c 4 = 




^) 






C5 = 


(c 2 - 


ci)detg, 


C6 = 


4 





•5+-fl r -+> 



We see that the critical temperature depends on non- 
magnetic disorder, but in contrast to the unconventional 
case, the effect is not described by a universal Abrikosov- 
Gor'kov function. At weak disorder the suppression is 
linear in the scattering rate, but with a non-universal 
slope: 



T r = T, 



c0 



ci 



C3 



C5 



7T 

8r' 



(77) 



In the opposite limit of strong impurity scattering, 
tT c0 < 1 (but still re c > 1, see Ref. H3), we use 
¥(x) = ln(l/x) + 0(1) > 1 at x -> 0, to find that the 
critical temperature approaches 



\9 2c 3 



(78) 



i.e. superconductivity is not completely destroyed by im- 
purities. The explanation is the same as in the conven- 
tional two-gap superconductors, see e.g. Refs. I49ll50ll5ll : 
Interband impurity scattering tends to reduce the differ- 
ence between the gap magnitudes in the two bands, which 
costs energy and thus suppresses T C1 but only until both 
gaps become equal. We have checked numerically that, 



for any positive-definite matrix g and for any p\ satisfy- 
ing the constraint p+ + p_ = 2, both the coefficient in 
front of t -1 in Eq. (|77[l and the exponent in Eq. ((78|) 
are negative, i.e. T* < T c q. 

In the model (|33[) . in which the pairing is isotropic 
and the order parameter has only one component, there 
is an analog of Anderson's theorem: The nonmagnetic 
disorder has no effect on the critical temperature, since 
the right-hand side of Eq. ([76]) is identically zero. The 
same is also true in the single-band limit, i.e. when the 
only nonzero constants are g ++ = g and p + = 2. 



B. Upper critical field 

To illustrate the effect of disorder on the upper critical 
field, we consider a tetragonal crystal, G = C^, in the 
field B — Bz, in which case H c2 (T) can be calculated 
analytically, see also Ref. [52|. The GL free energy 
takes the form 



T 



3 r^^(r){A 

w 

+ [k£ x ,(D 2 x + D 2 y ) + K\ y Dl] } V x(r), (79) 



where A, K^^K^ are real symmetric matrices, see Eqs. 
(f6"Tj) and (|68p in the conventional pairing case, and Eqs. 
(|65p and (|66p in the unconventional pairing case. 

In order to find the spectrum of the matrix differential 
operator <D\\< = K^ x , (D 2 X + D 2 y ) + k\ x ,D 2 7 , we introduce 
the operators 



he D x ± iD y 



(80) 



It is easy to check that a+ = a[_ and [a_,a+] = 1, and 
therefore a± have the meaning of the raising and lowering 
operators, while a 3 = a \ commutes with both of them: 
[ a 3i o±] — 0. Representing Ow in terms of the operators 
([80]). we have 

Oxx' - ^[K&,(4a+a- + 2) + K^a 2 }. (81) 

To calculate the eigenvalues explicitly, it is convenient to 
use the basis of the Landau levels \N,p), which satisfy 



a+\N,p) = VN+ l\N+ l,p) 
a-\N,p) = VN\N- l,p) 
a 3 \N,p)=p\N,p), 

where N — 0,1,..., and p determines the modulation 
along the z-axis: p — k zx /hc/eB. Writing the order 
parameter as a linear combination of the Landau levels, 



N,p 



r)x,N, P (r\N,p), 
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we obtain: 



(82) 



N,p XX' 



where 



Cxx'(N,p) 



K x L x ,(4N + 2) + Kl x ,p 2 



eB 

he 



(83) 



The upper critical field H C 2(T) can be found from the 
equation det C = 0, after maximization with respect to 
N and p. We assume that the maximum is achieved at 
N = p = and find the following expression for the slope 
of H c2 (T) at B -> 0: 



R = 



dH, 



<-2 



dT 
$ 



(84) 



T=T C 



A. 



- 2A + _a H 



2tt Af+Jf: 1 



A—K^ + - 2A+-K^_ 



T=T a 



where $o = ^hc/ e is the magnetic flux quantum, and 
Px&xx' - (px$XX' - \pxPX'{<l>x){4>X')} 

r c V: 



axx> 



2T~ r 



1 



1 



4tttT c 



4tttT c 



While the expression (f8"4")l is valid for arbitrary disor- 
der strength, we are especially interested in the limit of 
weak disorder. At tT c q 3> 1 the critical temperature can 
be represented as T c = T c0 [l — 6(7r/8rT c0 )], where the 
dimensionless coefficients b takes different values for con- 
ventional and unconventional pairing, see Eqs. (|77[) and 
(TT4")) . respectively. The matrices A, K^, a can also be 
expanded in powers of r _1 : 

A = A + SA, = + 5K L , a = a a + 5a, (85) 

where 

N F ( 1 



A 



XX' 



2 \i\p 
7ft 

0,XX' = ^A 



7cm 2 N F 2 2 

An >v = OAV 32^2^2 P^\ < rX v X,x) ^ 

N F 

a 0,AA' = nm PxSxX' , 



2T, 



c0 



is the Riemann zeta-function, and 



<L4 



AA' 



IT 



(1 - b)p\S\y - -pxPX'{<Px)(<l>X>) 



8rT c0 ' 



^AV 



42C(3)6\ 2 2 

-3H ^ I PX<>XX>{(p\V XtX ) 



+ 2/°APA'((^a)(0A'Va',x) + (</>A')(<Ma,z>) 



A A 



-TTPXPX' X! (li4 ) X'U\,x) {yifiX'VX' 



8tT c 



rO 



8 a 



AA' 



AA' 



[to get this we used C(4) = tt 4 /90]. 

In the clean case, Ao,aa' can be expressed in terms of 
the matrix elements of g, see Eq. (|69|) , which yields 



Rn — 



87T$nT f 



o^co 9+++9- 



2g- 1 detg 



7C(3)h 2 Exidxx-g^detg^lvl)' 



(86) 



with g defined by Eq. (|72[) . In the presence of impurities, 
we substitute the expansions ([55)) in Eq. (f8~4")) and obtain 
the correction to the upper critical field slope: 



SR 
R^> 



9++ +fJ — - 2.9 1 det g 

V2 

'ExifJxx-g-'detg^lvlJ 



8tT, 



rO 



(87) 



where 
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Vn 



\{p+9+- + p-g-+){<t>+){<t>-) - 

A=± 

7T 2 fy-r ( det 5^ |742£(3)& 
103)161 



9\\ ~ ( 1 + - ) det g 



( 



■A=± 



ffAA 



5 



V 7T 2 



^a,*) + Px((f>x)(<l>xvt,x) + y XI (t^a^a,*} 



x,y.z 



- det g 



2— x,y.z 



r 



It does not seem to be possible to draw any conclusions 
from Eq. ([57)1 about the effect of impurities on the slope 
of H C 2 in the general case, i.e. for arbitrary values of the 
coupling constants g\\> and the densities of states p\. 
For this reason, we just look at two limiting cases. 
In the model (|33|) . the expression ([86| takes the form 



Rn — 



167T$oT c0 



1 



(88) 



Neglecting the differences between the Fermi velocities 
and the densities of states in the two bands (these differ- 
ences are of the order of \j\/cf) and assuming a spher- 
ical Fermi surface, we obtain: p+{v\ x ) + p-{v^_ x ) — > 
2vp/3. In this way we recover Gor'kov's expression for 
the slope of the upper critical field in a clean isotropic 
superconductor^ In the same approximation, Eq. (|57|) 
becomes 



5R 



1 



168C(3) tT c{ 



(89) 



which coincides with the one obtained in Ref. [53| for the 
impurity-induced enhancement of the upper critical field 
in a conventional isotropic BCS superconductor. 

In the single-band limit, setting g\y = gS\, + 5\^ + , 
px = 2<5a i+ , v + (k) = v(k), </>+(fe) = (j>(k), we obtain 
from Eq. (|86|) : 



8^oT e0 1 
7C(3)ft2 (^2) 



and from Eq. (|H7|): 



8R .3 



1 



R 168C(3)rT c0 ' 
The coefficient C is given by 

r _ r _ , _ Si(7i^) 2 

in the conventional pairing case, and 

42C(3) Ei<7i<M< 



C — C lt 



(90) 



(91) 



(92) 



(93) 



in the unconventional pairing case. Since 42£(3)/7r 2 ~ 
5.115, C uncon is always negative, i.e. the upper crit- 
ical field for unconventional pairing is reduced in the 
presence of impurities. In the conventional case, the 
sign and the magnitude of the correction depend on 
the band structure. For example, consider a cylindri- 
cal Fermi surface with v(k) = VF(k x ,k y ,0)/k±, and 
■y(k) = (k y , ~k x ,0)/k± (the Rashba model), where k± = 

y/fc 2 + ky . Calculating the Fermi-surface averages in Eq. 

([92]) . we find C con = 1/2 > 0, which corresponds to the 
upper critical field enhancement by disorder. 



C. Application to CePt 3 Si 

Application of our results to real noncentrosymmetric 
materials is complicated by the lack of a definite informa- 
tion about the superconducting gap symmetry and the 
distribution of the pairing strength between the bands. 
One can only make progress by using some simple mod- 
els. For example, there are experimental indications that 
most of the Fermi surface in CePtaSi remains normal. 1 
It is plausible that the superconductivity resides in one 
of the two quasi-two-dimensional 7 bandsJ^ Therefore 
we take into account only the 7 bands and assume that, 
while the "+"-band is superconducting, the "-"-band re- 
mains normal, i.e. V-\ = V. =0. In contrast to the 

purely single-band limit, in the case under consideration 
P-^0- 

As far as the pairing symmetry is concerned, there 
is strong experimental evidence that the supercon- 
ducting order parameter in CePtsSi has lines of gap 
nodes i 54 i 55 ' 56 i 57 ' 58 The lines of nodes are required by sym- 
metry for all nontrivial one-dimensional representations 
of Cn v (A2, B\, and B2), so that the superconductiv- 
ity in CePtaSi is most likely unconventional. This can 
be verified using the measurements of the dependence of 
T c on the impurity concentration: For all types of un- 
conventional pairing, the suppression of the critical tem- 
perature is described by the universal Abrikosov-Gor'kov 
function, see Eq. (fT5|) . For the upper critical field, one 
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obtains from Eq. ((87)) : 

SR _ 
Ro 



C 



168C(3) rT c0 ' 



where 



C = 3 



42C(3) p+E ( (7i0+«+,*) 



2(0> 2 



(94) 



(95) 



One can see that for all types of unconventional pairing 
C < 0, i.e. the slope of H c2 is reduced by disorder. 

It should be mentioned that the lines of gap nodes can 
exist also for conventional pairing (A\ representation), 
in which case they are purely accidental. While the ac- 
cidental nodes would be consistent with the power-law 
behavior of physical properties observed experimentally, 
the impurity effect on T c in this case is qualitatively dif- 
ferent from the unconventional case. Indeed, using Eq. 
(|70| we obtain the following equation for the critical tem- 
perature: 



In- 



~Tc. 



from which we have 



T c = Too 1 - b 



(l-^ + } 2 )F(TT c ), (96) 



b=l-^ + ) 2 , (97) 



8rT d 



at weak disorder (tT c q 3> 1), and, using the asymptotical 
expression (f75|) . 



T r = T, 



<::0 



/ 7TrT ( 



cO 



V 7 



1. 



(98) 



at strong disorder (tT c0 -C 1). From the Schwarz in- 
equality (</>+) 2 < {(j>\) = 1, we find b > P-/2 and 
a > p-/p + . This means that anisotropy of the con- 
ventional order parameter increases the rate at which T c 
is suppressed by impurities. Unlike the unconventional 
case, however, the superconductivity is never completely 
destroyed, even at strong disorder. 



CONCLUSIONS 



is a mixture of spin-singlet and spin-triplet states, even 
without a spin-triplet term in the pairing interaction, see 
Ref. El). 

We have also studied the impurity effect on the criti- 
cal temperature T c and the upper critical field H C 2- Al- 
though T c is generally suppressed by impurities, this hap- 
pens differently for conventional and unconventional pair- 
ing: Any deviation of T c (r) from the Abrikosov-Gor'kov 
curve, see Eq. (|T3[) . in particular an incomplete sup- 
pression of superconductivity by strong disorder, is a sig- 
nature of conventional pairing symmetry. The impurity 
effect on the slope of H C 2 turns out to be sensitive to the 
pairing symmetry and the band structure, and therefore 
non-universal. 

In the general case, i.e. for arbitrary values of the intra- 
band and interband coupling constants and the densities 
of states, the microscopic expressions for the coefficients 
in the GL functional are rather cumbersome and there- 
fore of limited utility. Considerable simplifications occur 
only in some cases, for instance, in the BCS-like model 
(|33[) . whose properties resemble a conventional isotropic 
superconductor: The order parameter has only one com- 
ponent, T c is not affected by disorder, while H C 2 increases 
with disorder. On the other hand, in CePtaSi, where 
the pairing is likely unconventional, with lines of the gap 
nodes, we obtain that both T c and H C 2 are suppressed 
by disorder. 

As far as experiment is concerned, the upper criti- 
cal field measurements in CePtsSi, both in polycrystals 1 
and in single crystals^ give the zero-temperature values 
of H C 2 exceeding the Clogston-Chandrasekhar expression 
for the paramagnetic limit, Hp ~ 1.24T c //ib. This can 
be understood as being due to the l arge r esidual spin 
paramagnetic susceptibility, see Refs. Il4ll2lll22l . Surpris- 
ingly, the values of both the critical temperature and 
the upper critical field in polycrystals are higher than in 
single crystals. This is opposite to what has been ob- 
served in other unconventional superconductors and also 
disagrees with our theoretical predictions (assuming that 
the polycrystalline samples are intrinsically "more disor- 
dered" than the single crystals). To resolve this puzzle 
more work needs to be done on the purity dependence of 
the properties of noncentrosymmetric superconductors. 



We have derived the Ginzburg-Landau free energy of a 
noncentrosymmetric superconductor with isotropic non- 
magnetic impurities, using the microscopic model with 
a large spin-orbit splitting of the electron bands. If the 
pairing corresponds to a one-dimensional representation 
of the crystalline point group, then the order parame- 
ter has two components, one for each band, and the GL 
functional has the same form as for a two-band supercon- 
ductor (in the spin representation, the order parameter 
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